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Abstract 

We describe an extension of the nonlinear integral equation (NLIE) method to 
Virasoro minimal models perturbed by the relevant operator ^(1,3)- Along the way, 
we also complete our previous studies of the finite volume spectrum of sine-Gordon 
theory by considering the attractive regime and more specifically, breather states. 
For the minimal models, we examine the states with zero topological charge in detail, 
and give numerical comparison to TBA and TCS results. We think that the evidence 
presented strongly supports the validity of the NLIE description of perturbed minimal 
models. 
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1 Introduction 



Finite size effects play an essential role in the investigation of (both integrable and non- 
integrable) 1 + 1 dimensional quantum field theories QFT. They provide a possibility to 
determine many important physical characteristics such as S-matrices, mass ratios, re- 
lations between parameters appearing in the ultraviolet and infrared descriptions of the 
theory and a great deal of qualitative information on the spectrum. Besides their appli- 
cations in the study of low-dimensional QFT, they appear naturally in statistical physics 
and in the context of lattice field theory as well. 

In this paper we will focus on integrable theories. One of the first approaches to 
compute finite size effects in integrable theories was the Thermodynamic Bethe Ansatz 
(TBA) [|TJ which was used to calculate the vacuum (Casimir) energy 0]. The method was 
later extended to include ground states of charged sectors 0. More recently, using analytic 
properties of the TBA equations extended for complex values of the volume parameter, 
an approach to get excited states was proposed in 0|. Their method to get excited states 
sheds light on the analytic structure of the dependence of scaling functions on the spatial 
volume and up to now was the only method developed to deal with excited states in 
perturbations of minimal models. Its main drawbacks are that (1) it can be used only 
for systems for which a TBA equation describing the vacuum in finite volume is known; 
(2) to obtain the equation for a given excited state one has to do analytic continuation 
for each case separately, and a major part of this continuation can only be carried out 
numerically. Because of the requirement of the knowledge of the vacuum TBA equation and 
the complications of the analytic continuation, this method is limited at present to simple 
cases of integrable perturbations of Virasoro minimal models and some other perturbed 
conformal field theories. Similar results were obtained in || ||. 

This paper reports on a novel approach to the excited states of RCFTs in finite vol- 
ume, based on the nonlinear integral equation (NLIE) method, which has its origin in the 
so-called light-cone lattice Bethe Ansatz approach to regularize integrable QFTs. It was 
argued in [^] that sine-Gordon theory can be regularized using an inhomogeneous 6-vertex 
model (or equivalently, an inhomogeneous XX Z chain). The NLIE was originally devel- 
oped in this framework to describe the ground state scaling function (Casimir energy) in 
sine-Gordon theory in [||, || and it was shown that in the ultraviolet limit it reproduces 
the correct value of the central charge c = 1. We remark that similar methods were 



independently introduced in Condensed Matter Physics by other authors |10 



The NLIE was first extended to excited states in ||IT] where the spectrum of states 
containing only solitons (and no antisolitons/breathers) has been described. Using an idea 
by Zamolodchikov [|12|] they also showed that a twisted version of the equation was able 
to describe ground states of unitary Virasoro minimal models perturbed by the operator 
$(i j3 ). A framework for generic excited states of even topological charge in sine-Gordon 
theory was outlined by Destri and de Vega in (TH|. However, there has been a contradiction 



between the results of the two papers, which was resolved in |14|, [L5|] where we showed that 
it was related to the locality and the operator content of limiting ultraviolet conformal field 
theory (CFT). Besides that, we gave strong evidence for the correctness of the predicted 
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spectrum by comparing it to predictions coming from the truncated conformal space (TCS) 
method, pioneered by Yurov and Zamolodchikov in and extended to c = 1 theories 
in [14], Ol. Later we conjectured a modification of the NLIE to describe the states of 



sine-Gordon/massive Thirring theory with odd topological charge [|17 



The NLIE for sine-Gordon theory was generalized to models built on general simply- 
laced algebras of ADE type in |18| for the case of the vacuum. More recently, in |ll| P. 
Zinn- Justin extended the method to the spectrum of excited states for these models and 
he also made a first attempt to describe perturbations of minimal models of CFT. We will 
return to discussing his results and their relations to the present paper in section 4. 

The purpose of the present paper is to present a framework for describing general 
excited states of Virasoro minimal models (including nonunitary ones) perturbed by the 
operator $(1,3) (we consider only the massive case). In order to do that first we extend our 
previous studies to describe breather states in section 2. We show that in the IR limit the 
scaling functions resulting from the NLIE match the S-matrices conjectured in |20|] and 
then proceed to a numerical comparison to TCS data. We also discuss some consequences 
coming from the evaluation of the conformal weights of the UV limiting states. 

In section 3 we give a description of ground states of the perturbed minimal models and 
check them against TCS and, where available, TBA data. The comparison with TBA is 
especially powerful, since the numbers coming from the TBA method are exact up to the 
numerical precision of the iterative solution (which can be as small as computing power 
allows). Then in section 4 we write down and discuss the NLIE for the excited states. In 
our exposition we restrict ourselves to neutral states. In section 5 we check the simplest 
examples of the resulting excited state energies both numerically and qualitatively against 
TCS, and in section 6 we give our conclusions. 



2 The NLIE for the sine-Gordon theory in the attractive 
regime 

2.1 Notations and conventions 

Our convention for the sine-Gordon Lagrangian is 

C sG = j + A : cos (/?$) :J dx , 



(2.1) 



where (3 is the coupling constant and the dimensionful parameter A essentially defines a 
mass scale, which can be expressed in terms of the soliton mass M Ell] as follows: 
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(2.2) 



For later convenience, we define a new parameter p by 

In the repulsive regime p > 1 while in the attractive p < 1. The free fermion point is at 
p = 1 and the fcth breather threshold is p = 1/k. For later reference we recall that the 
UV limit of the theory is a c = 1 CFT, whose spectrum of primary fields consists of vertex 
operators V n , m with the conformal weights 

A ± = -(^±-mR) 2 , (2.4) 



2 \R 2 

where R is the compactification radius of the c = 1 free boson, related to p by 

1 p 
2^2 ~ p+ i • 

Let us briefly recall the NLIE for sine-Gordon theory. We will use the notations and 
conventions of the paper which the reader is invited to consult for more details. We 
put the sine-Gordon model on a cylindrical spacetime, with the infinite time direction and 
compact spatial extension of length (volume) L. The NLIE is a complex nonlinear integral 
equation for the counting function Z(i9): 



/oo 
dxG{§ -x-irj) log (1 + (-\) 5 e iZ{ - x+ir >)) 
-oo 

/oo 
dxG{d -x + irj) log (1 + (_i)<5 e -^-»?)) 
"OO 



(2.5) 



where / = ML is the dimensionless volume parameter and the parameter 5 takes the values 
or 1. The kernel of the equation is given by 



i r+oo 

Gi-ff) = — dke 



sinh w (P-V k 
ikti 611111 2 



2 sinh cosh f ' 



The function gifi^j) is the so-called source term, composed of the contributions from the 
holes, special objects (roots/holes) and complex roots which we call sources and denote 
their positions by the general symbol {$j} = {hk , Uk , Qc , Wk] (h stands for holes, y for 
special objects and c (w) for close (wide) complex roots). The source term takes the general 
form 

N H N s M c M w 

k=l k=l k=l k=l 

where 

= 2tt / dxG(x) (2.6) 
Jo 
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and the second determination for any function /($) is defined by 



{ )n \ /(#) -/(#- ip7rsign(9f77M?)) , p<l. ' 1 J 

The source positions i?j are determined from the Bethe quantization conditions 

Zfa) = 2*1, , Ij 6 Z + , (2.8) 

where I, are the Bethe quantum numbers. A complex root i?j is called close if j^m^l < 
min(7r, irp) , otherwise it is wide. Complex roots always come in complex conjugate pairs, 

p + 1 

except for self-conjugate roots which satisfy Qmfij = ±n — - — . Special roots/holes are 

real positions jjj where the Bethe quantization condition ( |2.8|) is satisfied, but the counting 
function is decreasing, i.e. 

Z' ( Vj ) < 0. 

Normally, Z is a monotonous increasing function on the real axis. 

C is an integration constant, which is a multiple of 7r and can always be set to zero by 
an appropriate redefinition of the counting function Z and the source term g (see (15|]). The 
redefinition affects only self-conjugate roots and may change the Bethe quantum numbers 
Ij from integer (i.e. 5 = 1) to half-integer (i.e. 5 = 0) or vice versa. 

The NLIE (|2.5|) only determines Z{$) in the fundamental analyticity strip |£ymi9| < 
min(7r,7rp). Outside the strip one must use the second determination of the counting 
function given by 

/oo 
dxG(p -x- vrj) n \og (1 + (-i)V z(x+M?) ) 
-oo 

/oo 
dxG{d -x + irj) n log (l + (-i)^-^"^) 
-oo 

(2.9) 

where C/j is yet another integration constant, which can be set to zero, similarly to C, by 
redefining the source terms appropriately. In addition, the Bethe Ansatz is periodic with 
a period m{p + 1), and so a fundamental domain for Z can be chosen as 

- 7F(P+I) <^< 7r(p+1) . 



The two determinations can be shown to suffice to cover the fundamental domain, outside 
of which Z is determined by periodicity. The counting function takes real values on the 
real axis and on the boundary lines of the fundamental domain, the so-called self- conjugate 
lines. 

The topological charge Q of any state can be obtained from the counting equation 

N H - 2N S = 2S + M C + 26(p - l)M w , (2.10) 
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where Nh is the total number of holes (ordinary/special), N$ is the number of special 
roots/holes, Mq is the number of close complex roots and Mw is the number of wide roots. 
S is called the total spin in the lattice Bethe Ansatz and is related to the topological charge 
by Q = 2S. It is known that Nh — 2Ns gives the number of solitonic/antisolitonic particles 
of the state [T^|. In the repulsive regime, the complex roots describe the internal degrees of 
freedom (polarization states) of solitons. In the attractive regime, however, configurations 
consisting entirely of wide roots describe the breathers and it is clear from eqn. fl2.10|) that 
they do not contribute to the topological charge. 

The energy and the momentum of a state can be expressed as 

N H N s 

E — Eb u ik = M cosh hj — 2M cosh yj 

3=1 3=1 
M c M w 

-M^coshcj - M^2(coshwj)n ( 2-11 ) 

3=1 3=1 

/°° dr 
—23m [sinh(x + irj) log(l + (-l) s e iZ ^ x+i ^)] , 
■oo 2vr 

N H N s 

P = MV sinh hj - 2M sinh yj 

3=1 3=1 

M c M w 

-M^sinhcj - M ^(sinh Wj) n ( 2 - 12 ) 

3=1 3=1 



r°° dx 

-M / — 2^sm [cosh(x + irj) log(l 



-1 °e 



S iZ(x+irj) 



where the values of hj, Cj, yj and Wj are fixed by the quantization conditions (|2.8|) . (Note 



that the sign of the wide root contribution is different from that in |T5fl - in our previous 



paper there is a misprint). The bulk energy term takes the form [^1 



E bu i k = —M 2 L tan ^. (2.13) 

For p and odd integer, this is divergent and is renormalized to give a term depending 
logarithmically on L (for an explanation see e.g. section 6 of (15|). The energy levels in 
finite volume take the general form 

E(L) = . (2.14) 

c(L) is called a scaling function and has the limiting value 

c(0) =c-12(A + + A~) , 
where are the conformal weights of the corresponding state in the UV limiting CFT. 



6 



2.2 Infrared limit and breather ^-matrices 



In the infrared limit I — > oo the term Zsinh($) develops a large imaginary part in the 
first determination away from the real axis, forcing the close complex roots to fall into 
special configurations called arrays (we use the terminology of [13|]). An array is a set of 
complex roots in which the roots are placed at specific intervals in the imaginary direction 
and have the same real part. In the attractive regime /sinh(^)// is nonzero and so this is 
true for nonselfconjugate wide pairs as well (in the repulsive case wide roots do not have 
such driving force), while self-conjugate roots have a fixed imaginary part anyway. The 
deviation of the complex roots from their positions in the array decays exponentially with 
I (see the paper |15] for an example with two holes and one complex pair in the repulsive 
regime). For a more detailed discussion, see |l3j. 

For the rest of this subsection, whenever it is not explicitly stated, we restrict ourselves 
to the attractive regime p < 1. The possible arrays fall into two classes: 

1. Arrays of the first kind are the ones containing close roots, which describe the polar- 
ization states of solitons. 

There are two degenerate cases: odd degenerate arrays, which have a self-conjugate 
root at 

. .7T(p+l) 

A, = + 1 2 

and accompanying complex pairs at 

nil - (2k + l)p) 



k = 0, 



2p 



and even degenerate ones, which only contain complex pairs, at the positions 



d k = ■& ± i 



7r(l - 2kp) 



fc = 0,... 



2p 



These arrays always contain exactly one close pair. The odd degenerate arrays in the 
repulsive regime reduce to single self-conjugate roots and the even degenerate ones 
to a single close complex pair. 

The description of the nondegenerate cases is a bit more complicated and can be 
found in They always contain two close complex pairs and they are the attractive 
regime analogous of wide pairs in the repulsive regime, but we will not need them 
here. 

2. Arrays of the second kind describe breather degrees of freedom. The odd ones contain 
a self-conjugate root 

O = + *7r(p+l)/2 

and wide pairs as follows: 

9m^ = t?±» 7r(1 - (2 f ! b + 1)p) , k = 0,...,s, 
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where 

while the even ones only contain wide pairs 

.tt(1 -2kp) 



1 , 



^5 



mdb = ± i , k = 0, 



and s runs in the same range. They correspond to the (2s + l)-th breather B 2s +i 
and the (2s + 2)-th breather B 2s +2, respectively. 

As one can see, arrays of the second kind become degenerate ones of the first kind, if 
we analytically continue increasing p. The reason is that breathers are of course soliton- 
antisoliton bound states, while degenerate arrays of the first kind describe scattering states 
of a soliton and antisoliton, as we will see shortly. 

In the infrared limit one can drop all terms containing the integral of 

log(l + (_i)* e ±i*(*±*j)) 

because they exponentially decay with I. One can therefore compute the energy and 
momentum contribution of a array of the second kind corresponding to the breather B s . 
The energy-momentum contribution turns out to be 

2M sin — (cosh tf,sinhtf) , (2.15) 

where 1? is the common real part of the roots composing the array. This is just the 
contribution of a breather B s moving with rapidity 1?. Arrays of the first kind do not 
contribute to the energy-momentum in the infrared limit, which lends support to their 
interpretation as polarization states of solitons. 

So far we have just reviewed some fundamental facts already known in the Bethe Ansatz 



literature (see flDJ and references therein). Now we proceed to show that with the above 
interpretation the NLIE correctly reproduces the two-body scattering matrices of sine- 
Gordon theory including breathers, which has not been done before. For the repulsive 



regime the S'-matrices were calculated in detail in our previous paper |T5| . 

Let us start with breather-soliton matrices. The Bethe quantization conditions for a 
state containing a soliton (i.e. a hole) with rapidity #1 and a breather B s with rapidity $2 
take the following form in the infrared limit. For the hole we get 

s 

Z{d x ) = Msinh^ - J2x($i - #2 - m)n = 27r/i, 

k=0 

where we denoted the prescribed imaginary parts of the roots in the array B s by pk- Here 
we used x(0) = to eliminate the source term for the hole. Now we can compute the 
second determination of x to be 
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gd($ + in/2) + gd{d - in/2 + in(p + 1)) , Smi? < -np 
gd{d - in/2) + gd{d + in/2 - in(p + 1)) , 3mi9 > np 
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where 

gd(0) = Z log . ry-r , 

sinn(wr/4 — v/ 2) 

which is essentially the Gudermannian arctan(sinh($)) with a suitable choice of branches 
T3|| . Now it is a matter of elementary algebra to arrive at 



Z(# x ) = M sinh tfx - i log 5 5Bs - tf 2 ) = 2tt/i 

where SsB a {$i — $2) is the soliton-breather S'-matrix conjectured in |2( 
One can start with the breather quantization conditions, too. Writing 

Z (tf 2 + ip k )n = Msinh (t? 2 + ipk) n +x(^2-A+ipk)ii+- ■ ■ = ^if ,k — 0,...,s (2.16) 

(the dots are terms due to wide root sources themselves, which cancel out up to multiples 
of 2ix in the next step) . Summing up these equations one arrives at 

2M sin ^ sinh(tf 2 ) - i log S SBa (0 2 - #1) = 2tt/ 2 , 

where 7 2 is essentially minus the sum of the quantum numbers of the wide roots composing 
the array (shifted by some integer coming from summing up the terms omitted in eqn. 

O). 

Using a similar line of argument we also reproduced the breather-breather S-matrices 
by writing down the Bethe quantization conditions for a state with two degenerate strings 
B s and B r of the second kind. One has to be careful that when Z("&) contains wide root 
sources which are expressed in terms of x(^)iii the second determinations of these terms 
will appear in Z(-d) n , i.e. terms that can be written roughly like (x($)ii)h ■ 

We remark that breather S-matrices of the sine-Gordon model were also reproduced in 
22| from the XX Z spin chain. However, their calculation is made on a homogeneous lattice 



XX Z chain, while we perform the derivation in the continuum limit of an inhomogeneous 
chain, which is described by the NLIE (|2.5|). It is very important that we can identify the 
spectral parameter $ with the rapidity of a relativistic breather particle with the correct 
mass, using eqn. ( |2.15| ), while in the case of the homogeneous lattice XX Z chain the 



dispersion relation is not the relativistic one. The relativistic dispersion relation ( |2.15|) can 



only be obtained by taking a special continuum limit, in which we send the inhomogeneity 
parameter to 00 and at the same time take the lattice spacing to (for details see |T3| ). 

Scattering state of a soliton and an antisoliton can be described by taking two holes 
and a degenerate array of the first kind. There are two possibilities now, corresponding to 
scattering in the parity-odd and parity-even channels. Following the procedure outlined 
in fl~5|| , we were once again able to reproduce the corresponding scattering amplitudes. 
The results presented here together with those of exhaust all two-particle scattering 
amplitudes of sine-Gordon theory, both in the repulsive and attractive regime. 
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/ 


TCS 


NLIE 


TCS 


NLIE 


1.0 


2.38459 


n/a 


1.84996 


n/a 


1.5 


1.68168 


n/a 


1.30438 


n/a 


2.0 


1.35391 


n/a 


1.05038 


n/a 


2.5 


1.17420 


n/a 


0.91152 


n/a 


3.0 


1.06692 


1.0655810539 


0.82903 


0.8285879853 


3.5 


0.99985 


0.9980153379 


0.77773 


0.7771857432 


4.0 


0.95664 


0.9542867454 


0.74499 


0.7443106400 


4.5 


0.92845 


0.9254766107 


0.72381 


0.7229895177 


5.0 


0.91000 


0.9063029141 


0.71006 


0.7090838262 



Table 1: The first breather state at p = — and p = —. Energies and distances are measured in 

7 9 

units of the soliton mass M, and we have subtracted the predicted bulk energy term from the 
TCS data. 



2.3 Some examples of breather states 

The vacuum scaling function and the multi-soliton states of sine-Gordon theory (both in 
the attractive and the repulsive regime) have already been examined in jrj|, where we 
found agreement with TCS predictions. So we proceed to take a look at the simplest 
neutral excited state, which is the one containing a first breather B 1 at rest. The source 
term to be written into the NLIE turns out to be 

cos — — i sinh $ 
g($) = -i log ± , (2.17) 

cos h % sinh ■d 

2 

and the self-conjugate root is located exactly at t?o = wr(p + l)/2, since the breather has 
zero momentum. There is no need to look at the Bethe quantization condition as the root 
does not move due to the left-right symmetry of the problem. This state is quantized 
with integer Bethe quantum numbers, i.e. 5=1 (cf. fli~5l |l7]] for an explanation of the 
connection between the choice of 5 and locality of the corresponding quantum field theory). 



Calculating the UV conformal dimensions along the lines presented in |T5| we obtain 
the conformal dimensions 

A± = p+T ' < 2 ' 18 > 

so the ultraviolet limit of this state is a linear combination of the vertex operators V±i o of 



the c = 1 UV CFT (see eqn. (|2.4j) ). This is in perfect agreement with the TCS calculations 
performed by us (for an outline of the method see |15|). To normalize the TCS, we used 
the formula ( |2.2|) . 

Table [j] presents the energy values obtained by iterating the NLIE in comparison to 
results coming from TCS, at the values p = 2/7 and p = 2/9. 

The table shows that iteration of the NLIE fails for values of I less than 3 (the actual 
limiting value is around 2.5). What is the reason? 
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Figure 1: The behavior of the function Z{p) at / = 3 and I = 2, respectively. 



We plot the counting function Z on the real line for I = 3 and I — 2 in figure [p. In a 
first approximation we can safely neglect the integral term for these values of the volume to 
see the qualitative features that we are interested in. What we see is that the behaviour of 
the function changes: its derivative changes sign at the origin. As a result, two new holes 
appear where the new real zeros of the function are. But the topological charge remains 
zero, due to the fact that now we have a special root at the origin and so N s = 1 and 
Nh = 2. The two new holes do not give us any new dynamical degrees of freedom: their 
quantum numbers are fixed to be and so their positions are uniquely determined. 

Of course when we calculated the UV dimension, the appearance of the new sources 

had to be taken into account to obtain the result ( |2.18| ). It turns out that for - < p < 1 the 

two holes are left/right movers, while for p < — they remain central, using the terminology 

3 



of 15 



How does this phenomenon affect the iteration scheme for the NLIE? The two new 
zeros of which is a complex analytic function apart from logarithmic branch cuts, 

actually correspond to singularities of the logarithmic term in the NLIE (|2.5|) . They come 
along the imaginary axis in the $ plane as we decrease /, and at a certain point they cross 
our integration contour which runs parallel to the real axis at distance r/. As they make 
the logarithmic term in our NLIE ( [2.5| ) singular, they blow up our iteration scheme. After 
reaching the origin of the plane (at exactly the radius where the derivative of Z becomes 
0), they continue to move along the real axis, which corresponds to crossing a square root 
branch cut. 

We do not go into details here as this problem is currently under investigation^. We just 
remark that these issues prove to be highly nontrivial and for the time being, unfortunately, 
they prevent us from having a reliable numeric scheme for the NLIE below the critical 
volume. The problem will be examined in a forthcoming publication. 

One can estimate the volume where the slope of the counting function changes sign by 

4 Work in progress in collaboration with P. E. Dorey and C. Dunning, Durham. 
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neglecting the integral of the logarithmic term. The result is 



I critical nryfl (^■^^) 

cos T 



which gives a value of around 2.22 for p = 2/7 and 2.13 for p = 2/9. The actual limiting 
value is a bit higher, partly due to the finite value of rj used in the iteration program and 
partly because the iteration already destabilizes when the singularities come close enough 
to the contours. It must also be noted that the integral term cannot eventually be neglected 



when the singularities are close to the contour, which is an additional reason why (|2.19|) is 
just a crude estimate. 

We make a short digression to examine the UV limit of the second breather Bi. The 
second breather at rest is described by a wide pair at positions 

i? = ±i- . 
2 

Calculating the UV conformal dimension we get 

P 



A+ = A" 



P+ 1 



which turns out to be the same as that of the first one ( |2.18|) , i.e. this state must orig- 
inate from the other linearly independent combination of the vertex operators V±\$ in 
the ultraviolet. This is again in perfect agreement with TCS and confirms a result by 
Pallua and Prester [^] who used XXZ chain in transverse magnetic field to regularize 
sine-Gordon theory. They calculated scaling functions numerically on a finite lattice for 
several concrete values of p, and arrived at this conclusion by looking at the numerical 
data. However, our method to compute UV dimensions gives us an exact analytic formula 
and therefore much stronger evidence. This result is interesting because it invalidates a 



conjecture made previously by Klassen and Melzer |]24[ who identified the second breather 
as a linear combination of V±2,o- 

To close this section, we present the lowest lying example of a two-breather state, 
containing two Bi particles with zero total momentum. It turns out that this is a state for 
which the numerical iteration of the NLIE is not plagued with the problem found above 
for the first breather. Locality constrains the state to be quantized with half-integers and 
for lowest energy the quantum numbers of the self-conjugate roots must take the values 

/ 1 = I, a I 

1 2 ' 2 

We remark that in contrast to the case of holes, the self-conjugate root with / > moves 
to the left, while the one with / < moves to the right. This is due to the fact that 
the second determination of Z is in general a monotonically decreasing function on the 
self-conjugate line. In order to determine the position of the two self-conjugate roots we 
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/ 


TCS 


NLIE 


1.0 


12.1601 


12.159257 


1.5 


8.20139 


8.2006130 


2.0 


6.24771 


6.2465898 


2.5 


5.09489 


5.0937037 


3.0 


4.34132 


4.3397021 


3.5 


3.81513 


3.8129798 


4.0 


3.43020 


3.4275967 


4.5 


3.13912 


3.1357441 


5.0 


2.91308 


2.9089439 



Table 2: The two-breather state at p = =. 



need the second determination of Z. The second determination of the self-conjugate root 
source turns out to be 

i sin up — sinh i? 

ii = * log — . r q • 

z sin 7rp + sinh 17 

Up to some signs, this is just the phase shift which arises when two breathers scatter on 
each other, which is exactly why the IR analysis gives the correct scattering amplitude. 
Using this formula, we obtained the numerical data presented in table ||. The UV limit 
for this state can be calculated from NLIE to be a symmetric first level descendent of the 
vacuum with weights 

A + = A~ = 1 , 

which agrees with TCS. (Note that this descendent exists due to the fact that there is 
a U(1)l x U(1)r Kac-Moody symmetry at c = 1: this state exactly corresponds to the 
combination of the left and right moving currents J J.) 

We remark that the precision of c = 1 TCS is at maximum 5 — 6 digits for low lying 
states and small values of I, and the truncation error increases with the volume and as 
we go higher and higher in the spectrum. We also mention that to achieve this precision 
it is necessary to go up to around 4000 — 5000 states, which is a highly nontrivial task 
accomplished only by writing all the calculation in compiled C programs (including the 
evaluation of CFT matrix elements) and it was really stretching the computing power 
available to us to its limits. In contrast, for states without moving roots the NLIE can be 
readily iterated up to 12 digits precision by a simple personal computer and even for states 
with moving roots like the two-breather one it is not very hard to achieve a precision of 
6 — 7 digits (although the number of necessary iterations grows when one decreases the 
volume parameter I and the prescribed precision). 
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3 Zamolodchikov's a-twist and ground states for mini- 
mal models perturbed by $(13) 

Some time ago, Zamolodchikov has put forward the idea of modifying sine-Gordon theory 
by a twist a [|l^|. The NLIE describing this situation has the following form: 



/oo 
dxGi'd - x - w?) log (l + e lZ{x+iv) ) 
■00 

/oo 
dxG{d -x + ir])\og(l + e -«*G*-*»)) 
■00 



(3.20) 



For definiteness we choose half-integer quantization rule with 5 = 0, since it is obvious that 
shifting a by it and then redefining Z one can change the Bethe quantum numbers from 
half-integers to integers. The energy level determined by this equation has the leading UV 
behavior 



l-^rf-l ■ (3.21) 



with 



P + 1 \7T. 

Furthermore, it is well-known that the perturbation of the Virasoro minimal model Vir(r, s) 
by its relevant primary operator $(1,3) is integrable and is described by an RSOS restriction 
of sine-Gordon theory |25[ with 

p=— !— . (3.22) 

s — r 

We will use for this model the shorthand notation Vir(r, s) + $(1,3). 
Now, putting a = n/r we get 

c-=i-i, 

rs 

which is exactly the effective central charge of the minimal model Vir(r, s). Therefore one 
can expect that the twisted equation describes the ground state of the model Vir{r,s) + 
$(i,3)- In fact, Fioravanti et al. calculated these scaling functions for the unitary case 
s = r + 1 and showed that they match perfectly with the TBA predictions already available 
. Moreover, choosing the following values for the twist 

A/7T 

a = ± — , k = l...r- 1 (3.23) 
r 

they obtained the conformal weights of the operators , k = 1 ... r — 1 in the UV 

limit (the sign choice is just a matter of convention). In our notation, denotes the 

primary field with conformal weights 

A+ = A- = (*»-fr) 2 -(*-rf 

Asr v ' 
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The models Vir(r, s) + $(1,3) have exactly r — 1 ground states. In fact, one can see from the 
fusion rules that the matrix of the operator $(1,3) is block diagonal with exactly r — 1 blocks 
in the Hilbert space made up of states with the same left and right primary weights. In 
each of these blocks, there is exactly one ground state and for the unitary series s — r + 1, 
it was conjectured in |^6| that their UV limits are the states corresponding to <&(k,k)- One 
can check that in the general nonunitary case the twists ( |3.23Q correspond in the UV limit 



to the lowest dimension operators among each of the r — 1 different blocks of primaries 
(see explicit examples later). These ground states are degenerate in infinite volume, but 
for finite I they are split with their gaps decaying exponentially as / — > 00. In the unitary 
case, they were first analyzed in the context of the NLIE in [|ll]] where it was shown that the 
NLIE predictions perfectly match with the TBA results already available for the unitary 
series. 

However, ground states for nonunitary models have not been treated so far and therefore 
now we proceed to give examples of ground state calculations for nonunitary models. The 
models we select are the ones that will be used for comparison in the case of excited states 
as well. The first is for the scaling Lee- Yang model Vir(2,5) + $(1,3), for which we have 
also given data from TCS |16|] and TBA [|] for comparison (table ||). The notation used is 



l B = M B L , 

where 

M B = 2Msin — = V%M 

is the mass of the fundamental particle of the Lee- Yang model (this is more natural here 
than using the mass M of the soliton of the unrestricted sine-Gordon model as a scale, 
since the soliton disappears entirely from the spectrum after RSOS restriction). There is 
only one independent value of the twist, which we choose to be 

71 

a = — . 

2 

Here and in all other subsequent calculations the TCS data were normalized using the 



analogue of the coupling-mass gap relation (|2.2|) from 21]. There is only one ground state 



in this model, which corresponds to the primary field with conformal weights 

A+ = A- = -\ , 

5 

which is in agreement with TBA and TCS predictions. We have also found a perfect 
agreement for the models Vir(2, 7) + $(1,3) and Vir(2,9) + $(1,3), but we do not present 
those data here. We remark that the TCS for the minimal models converges much better 
than the one for c = 1 theories: all TCS data in table [3| and subsequent ones were produced 
by taking a few hundred states and in some fortunate cases (e.g. the ground state of the 
scaling Lee- Yang model for small values of I) we were able to produce data with up to 9— 10 
digits of accuracy! The better convergence meant that all the computation could be done 
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h 


TCS 


NLIE 


TBA 


0.1 


-2.0835015786 


-2.0835015787 


-2.0835015786 


0.5 


-0.3803475256 


-0.3803475281 


-0.3803475281 


1.0 


-0.1532068463 


-0.1532068801 


-0.1532068801 


1.5 


-0.0763483319 


-0.0763484842 


-0.0763484842 


2.0 


-0.0406269362 


-0.0406273676 


-0.0406273676 


2.5 


-0.0222292932 


-0.0222302407 


-0.0222302407 


3.0 


-0.0123492438 


-0.0123510173 


-0.0123510173 


3.5 


-0.0069309029 


-0.0069338817 


-0.0069338817 


4.0 


-0.0039198117 


-0.0039244430 


-0.0039244430 


5.0 


-0.0012721417 


-0.0012816882 


-0.0012816882 



Table 3: The vacuum of the Virasoro minimal model Vir(2, 5) perturbed by $(1,3) ■ The energy 
and the volume are normalized to the mass of the lowest excitation, which is the first breather of 
the unrestricted sine-Gordon model. The TCS data shown have the predicted bulk energy term 
subtracted. 



with the computer algebra program Mathematica, greatly simplifying the programming 
work. 

Models of the class Vir(2, In + 1) + $(1,3) have only one ground state. For models with 
two ground states, we can take a look at Vir(3, 5) and Vir(3,7). In the first case the 
ultraviolet spectrum is defined by the following Kac table, where the weight ( |3.24| ) of the 
field $(fc,i) is found in the k-th row and Z-th column. 






1 


1 


3 

4 


4 


s 


-V- 

20 






The two blocks of the perturbing operator $(1,3) are defined by the fields {$(1,2) , $(1,4)} 
and {$(1,1) , $(1,3)}, respectively. The ground states correspond in the UV to the operators 
$(1,2) and $(1,1) , as can be checked directly using formulae (|3.21|) , (|3.22|) and fl3.23j ). 

We also have TBA data to compare with, using the TBA equation written by Christe 
and Martins [27]]. The lower-lying ground state is obtained directly from their TBA, while 
for the other we used Fendley's idea of twisting the TBA equation [[]]. The numerical 
results are presented in tables | and ||. 

In the case Vir(3, 7) (table |^) we can only have a comparison with TCS results, but it 
still looks pretty convincing. The spectrum of the UV theory is 






5 


1 


3 


4 


5 










4 


5 













4 


7 


28 


7 


28 



Here the two blocks are {$(1,2) , $(1,4) , $(1,6)} and {$(1,1) , $(1,3) , $(1,5)}, respectively. 
The ground states correspond in the UV to the operators $(1,2) and $(1,3) . 

To summarize, we now have sufficient evidence to believe that the a-twisted NLIE 
describes the correct scaling functions for ground states of minimal models perturbed by 
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/ 


TCS 


NLIE 


TBA 


0.1 


-3.074916 


-3.0749130189 


-3.0749130190 


0.3 


-0.944161 


-0.9441276204 


-0.9441276204 


0.5 


-0.509764 


-0.5096602194 


-0.5096602194 


0.8 


-0.265436 


-0.2651431026 


-0.2651431026 


1.0 


-0.186038 


-0.1855606546 


-0.1855606546 


1.5 


-0.087300 


-0.0861426792 


-0.0861426792 


2.0 


-0.045910 


-0.0437473815 


-0.0437473815 


2.5 


-0.026746 


-0.0232421927 


-0.0232421927 


3.0 


-0.017868 


-0.0126823057 


-0.0126823057 


4.0 


-0.013546 


-0.0039607326 


-0.0039607326 



Table 4: One of the two ground states of the Virasoro minimal model Vir(3, 5) perturbed by 
$(i,3), corresponding to a = |. The energy and the volume are normalized to the mass of the 
lowest excitation, which is the soliton of the unrestricted sine-Gordon model. The TCS data 
shown have the predicted bulk energy term subtracted. 



/ 


TCS 


NLIE 


TBA 


0.1 


3.117844 


3.1178476855 


3.1178476853 


0.3 


0.985360 


0.9853990810 


0.9853990810 


0.5 


0.540427 


0.5405470784 


0.5405470784 


0.8 


0.282725 


0.2830552991 


0.2830552991 


1.0 


0.197143 


0.1976769278 


0.1976769278 


1.5 


0.089277 


0.0905539780 


0.0905539780 


2.0 


0.042960 


0.0453290013 


0.0453290013 


2.5 


0.019978 


0.0238075022 


0.0238075022 


3.0 


0.007209 


0.0128843786 


0.0128843786 


4.0 


-0.006592 


0.0039866371 


0.0039866371 



Table 5: The other ground state of the Virasoro minimal model Vir(3,5) perturbed by $(1,3), 
corresponding to a = j. The energy and the volume are normalized to the mass of the lowest 
excitation, which is the soliton of the unrestricted sine-Gordon model. The TCS data shown have 
the predicted bulk energy term subtracted. 
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ground state 1 


ground state 2 


1 


TCS 


NLIE 


TCS 


NLIE 


0.1 


-3.7039095315 


-3.7039095318 


0.7770591073 


0.7770591068 


0.5 


-0.6377729020 


-0.6377729703 


0.2075541041 


0.2075539646 


1.0 


-0.2314222239 


-0.2314229899 


0.1136501615 


0.1136487561 


1.5 


-0.1038698195 


-0.1038727861 


0.0649219172 


0.0649170326 


2.0 


-0.0505812388 


-0.0505886086 


0.0368105709 


0.0367993127 


2.5 


-0.0258481751 


-0.0258625389 


0.0208408686 


0.0208198502 


3.0 


-0.0136578583 


-0.0136819143 


0.0118382379 


0.0118037311 


3.5 


-0.0073877811 


-0.0742408962 


0.0067682934 


0.0067164173 


4.0 


-0.0040552779 


-0.0410602666 


0.0039103368 


0.0038372749 


5.0 


-0.0012235029 


-0.0013070235 


0.0013931394 


0.0012676541 



Table 6: The two ground states of the Virasoro minimal model Vir(3, 7) perturbed by 3>(i,3). 
The energy and the volume are normalized to the mass of the kink, which is the soliton of the 
unrestricted sine-Gordon model. Ground state #1 is obtained by putting a = ^, while for ground 
state #2 a = The TCS data shown have the predicted bulk energy term subtracted. 



^(1,3) even in the nonunitary case. However, the NLIE for sine-Gordon is known to work 
for excited states as well. But how do we get the excited state spectrum of the minimal 
models now? 



4 Excited states in minimal models perturbed by $(1,3) 
4.1 The excited state equation 

In this section we write down the NLIE for excited states of Virasoro minimal models per- 
turbed by $(1,3) and we also point out some essential differences from the results obtained 
by P. Zinn- Justin in ||19|| . 

The a-twist can be recast in the language of the original light-cone lattice Bethe Ansatz 



as a twist angle uo in the Bethe Ansatz equations (BAE) which take the form [28 



^ sinh — 

7T 




+ 6 + 


in 
~2 _ 


sinh — 

7T 


*i 


-6 + 


in 

y_ 


\ 


sinh — 
V vr 




+ 0- 


in 
~2 _ 


7 

sinh — 

7T 




-e- 


in 
~2 _ 


J 



N 



M sinh - pi?,- - -d k + in] 

-2iu> 7T 



k= i sinh — \dj — i9 



in 



n 



(4.25) 



where 



7 



n 



p+l 



and for the rest of the notations we refer the reader to [13, 15 
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From the BAE ( fl.25| ) the NLIE for excited states can be derived by the well-known 
methods [13|, We obtain the result 



/oo 
dxG{d -x-irj) log (l + (-l) s e iZ( - x+i ^) 
-oo 

/oo 
dxG{d -x + irj) log (1 + (-l) s e - iZ ^-^) 
-oo 



The parameter a turns out to be 



a 



p+1 

U V Xoo 

V 



s 



p + l 



10 
7T 



1 

2 + p 



5 



7T 



(4.26) 



(4.27) 



where 



Xc 



7T J9 

2~ 



P 



and S is the Bethe spin as determined from eqn. (|2.10| ). The square brackets denote the 
integer part, i.e. [x] is the largest integer smaller or equal to x. The relation ( [4.27] ) can be 
obtained in the following way. We recall that it is the derivative of the NLIE that can be 
derived from the lattice Bethe Ansatz using a contour integral trick [T3|, [TI|. Therefore it 
has to be integrated to give the NLIE itself and so there appears an integration constant, 
which can be determined by matching the asymptotic values of the lattice counting function 



Z coming from the NLIE to the asymptotic values determined from the Bethe Ansatz QT3 
This is the same procedure that yields the integration constant C in (|2 



]). However, while 

C is a multiple of n and can be absorbed by redefining the counting function Z and the 
Bethe quantum numbers Ij in ([Op (and possibly changing 5 to 1 — 5), this is not true for 
a as it can take any real value in general. 

the form of a is 



In contrast, in IS 



UJ 



a 



1 - 



7 



7T 



u- 



p +1 

V 



(4.28) 



and the choice for u is 7. However, this formula misses the second term in eqn. (|4.27| ). 
This term is extremely important since it guarantees that the formula (f4.27j) satisfies the 
property 

a — ► a + 2tt when u — > uj + n . 



This relation is required for self-consistency: since the BAE (|4.25|) are invariant under the 
shift of uj by tt, the NLIE (which is equivalent to the BAE before taking the continuum 
limit) must be invariant under this shift as well. (Note that shifting a by 2n is an invariance 
of the NLIE (14.26 ), with an appropriate redefinition of Z and the Bethe quantum numbers 

h •) 

From now on we restrict ourselves to the case of neutral (i.e. S = 0) states. It is 
easy to see that even for states with a zero charge the relation between a and uj is highly 
nontrivial. 
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Now we proceed to show that choosing 

u = k'y 

where k is integer, as the value of uj we can reproduce all the required values of a listed in 
equation (|3.23j ). First of all, we substitute the value of p from (|3.22|) to obtain 



UJ 



k(s — r)n 



Since r and s are relative primes, the independent values of uj mod n can be written as 

In 

uj = — , I = 0, . . . , s — 1 . 

s 

For S = 0, we can rewrite the formula ( H-271 ) as follows 



In 2r — s 

a = — + — n 

r 2r 



1 i 

2 + S 



1 / 

2 ~ s 



We are interested only in the value of a mod it, since using the parameter 5 one can 
effectively shift a by n. This leaves us with the formula 



In sn 
r 2r 



1 I 
2 + S 



1 I 
2~~s 



The first possibility is that I < f , which simply gives us the values 



In 



a 



When s is even, we can have I — |, which gives us a = 0. Finally, when / > |, we get the 



values 



a 



(l-s)n 



It is easy to check that these formulae reproduce every value 



nn 

a = — mod n 

r 



at least once, using the fact that s > r and that the values above form an uninterrupted 
sequence of s numbers (or when s is even, of s — 1 numbers, the zero repeated) with equal 
distances -. 

r 

As we have already seen in the previous section, all the values 



kn 

a = — , k — 1, 
r 



, r — 1 
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are necessary to reproduce correctly the r — 1 ground states of the model Vir(r, s) + $(1,3). 
The single choice 



7T 

a = — 
V 



of the paper jl9| is not enough to reproduce all the possible states. The twisted lattice 
Bethe Ansatz was analyzed by de Vega and Giacomini in (28|. On the lattice, passing 
from the sine-Gordon model to the perturbed Virasoro model amounts to going from the 
six-vertex model to a lattice RSOS model. In |^] it was shown that to obtain all the states 
of the RSOS model it is necessary to take all the twists 

uj = k r ) mod it 

into account. The fact that not all these twists correspond to inequivalent values of a and 
so to different physical states is a consequence of the RSOS truncation. 

To close this section we remark that the parameter a drops out of the second determi- 
nation of Z ( |2.9|) in the attractive regime. This is important because as a consequence the 
IR asymptotics of the breather states does not depend on a and so the S'-matrices involv- 
ing breathers are unchanged. In fact, scattering amplitudes between solitons and breathers 
remain unchanged too, as can be seen from examining the argument that we used to derive 
them in section 2. This matches with the fact that the RSOS restriction from sine-Gordon 
theory to perturbed minimal models does not modify scattering amplitudes that involve 



two breathers or a breather and soliton 25 . 



4.2 The UV limit 

There is in fact a very simple intuitive argument to show that the states we get from 
( |4.26| ) are related to the minimal models in the UV limit. To keep the formulae simple, we 
present it here for the unitary case p — r — s — 1. Let us first recall that the UV limit of 
the original NLIE (|2.5| ) yields the vertex operators Vc n , m ) and their descendants, where for 
general choice of 5, n is a half-integer and m is an integer. 

Let us look at the weights in the neutral sector, which means m = 0. Using the formulas 



for the UV limit of the NLIE from our previous paper ||15|| , one can see that introducing 
a 7^ is equivalent to shifting the quantum number nton + ^. One has to be careful that 
since the value of the central charge is shifted from 1 to the one of the minimal model, we 
have to take this shift into account when computing the conformal weight from the leading 



UV behaviour of the energy level ( |2.14| ) . We put in the value 



* = l,...,r-l 

V 



and the resulting conformal weights take the form 



. 1 . (2np + l) 2 -l 

A + = A-= v / . J — . 4.29 
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Comparing to the formula (|3.24j) we see that this is the weight of the field $(;,;_2n) in 
the minimal model Vir(p,p + 1). Therefore one expects that we get something related to 
minimal models, however, in order not to overflow the Kac table, the range of n must be 
restricted as 

1 < I - 2n < p . 

For charged states, a similar calculation can be performed. We find in particular that 

a: 
7T 



2 (A+ - A - ) = m- mod 1 , (4.30) 



so these states have fractional Lorentz spin in general. 

In a more precise way, the UV limit of the twisted NLIE can be examined with exactly 
the same method as outlined in |13|, so here we just give the formulae necessary to 
do the computations and for the derivation and notational details we refer to the papers 
above. The right/left kink equation can be obtained by substituting 

2 

— ?9 ± log - 

into the NLIE and keeping only the leading terms when I — > 0. It takes the form 

/oo 
dxG{d - x)Q±{x) , (4.31) 
-oo 

where the source term g±($) is the limit 

g±(&) = gU±\og?\ , I - 0, (4.32) 

and 

I _|_ ( — \y e iZ{x+iV) 

Q±(x) = lim — 2 log- 



77^+0 1 + ( — iy e -iZ{x-i V ) 

We call a source right/left moving if its position given in this limit by 

2 



^±l0gy 



where $ is finite as / — > oo. We denote the number of right /left-moving holes by N H and 



their positions by ft, • and similarly introduce the notations Ng , for the special objects, 
Mq , cf for the close complex roots, and M^> , ti^for the wide complex roots. There can 
also be sources whose positions remain finite as / — > oo : they are called central. We 
introduce as 

S ± = \\ n h- 2^5 -M%- 2M±6{p - 1)] . 



The plateau equation takes the form: 



Z±{Too) = a + #±(=F°o) + — u± , X oo = t:- — -, (4.33) 

7T I p 
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where we have 

Z±{ z fOo) = 0J± + TT5 + 27tk± , — 7T < U)± < 7T , 

and 

g±(=foo) =±2(S- Xoo + 2vrfc± . 

The numbers k^ depend on the configuration of wide roots and are generally integers 
except when we have an odd number of self-conjugate roots, in which case they are half- 
integers. From these relations one can compute the value of ui± to obtain 

W± = -^L ± 2. (S - 2S±) P -^r + 4x-*-k± , 
p+1 p + 1 I + p 

where k± are some (half)-integers. The conformal weights themselves can be computed 
using the formula 



24 v H s c WJ 2ir p 16tt 2 



(4.34) 



where 1^ is the sum of the Bethe quantum numbers of left/right holes etc., c is the central 
charge 

6(r — s) 2 
c = 1 , 

rs 

and the quantities E± are given by 



N i M c M v 



S ± = -EM^)+ 2 E s± (yf) + E ^ ( c ?) + E ^ K) « * 2s 
j=i j=i j=i j=i 



which can be written more explicitly as follows 

S ± = T 2S ± a - AS ± (S - S^Xoo + 2>irq± , (4.35) 

where is an integer or half-integer which depends on the configuration of wide roots. 
One has to be careful that g±($)n is defined by 



g±{#) n = g[#±log~) ,1^0, (4.36) 



and not as the second determination of <7±($). 

We close this section with some immediate consequences of the above formulae. First 
note that because the value of a for a minimal model is never a multiple of ir, one does 
not expect central sources in the UV limit (in all the examples of (15| with central sources, 
the left-right symmetry of the NLIE was crucial. This symmetry, however, only holds for 
a = or 7r). As a consequence we have 

S —— S ~\~ s , 
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and in addition u + = which means that there are only one-plateau systems. By an 
inspection of the formula Q4.34Q this implies that for any state with S = 



A+ - A" 



is integer or half-integer. In fact, choosing the quantization rule and the parameter 5 in an 
appropriate way, one can ensure that this difference is integer (see |TJ|, [17]]). This means 
that the UV limit of any neutral state is either a field occurring in the ADE classification 
of modular invariants |29|] or (in case we choose 8 so that A + — A~ is half-integer) it is a 
field from a fermionic version of the minimal model |24 



For the case of neutral states, an elementary calculation gives the following formula 



A+ = A" 



c- 1 
~2A 



+ 



(^- + (2k-AS+)p 
4p(p + 1) 



+ N. 



± 



where k, N± are integer or half-integer which upon substituting the value 



a 



In 
P 



is just identical to the formula ( |4.29|) , apart from N±, which must in fact be nothing else 
but descendant numbers and so one expects that they are eventually integer. A similar 
result was obtained in ItlSfl , but just as we said above the choice made in that paper does 
not allow one to fill the Kac table: for unitary models it produces only the fields 3>(i, n )- 
However, similarly to the argument presented at the beginning of this section, it is not 
clear whether the weights actually stay inside the Kac table, for which in the unitary case 
one must require 

1 < / - 2k + AS + < p . 

Due to the fact that the configuration of sources in the UV may be very non-trivially 
related to the one in the IR, this condition is very hard to check in general, but no concrete 
examples that we calculated have ever violated this bound. 

From the formulae (|4.34|) and ( [4. 35| ) it is also clear that in general 



2 (A+ - A - ) = 2S- mod 1 



and so we see again that general charged states will have fractional Lorentz spin. Remem- 
bering the relation m = 2S ||15j| . this agrees with our previous result ( |4.30|) . Actually, it is 
known that charged states in the models Vir(r, s) + $(1,3) generally have fractional Lorentz 
spin pO 



5 Concrete examples of excited states 



24 



h 


TCS 


NLIE 


0.1 


23.05277 


n/a 


0.5 


4.679779 


n/a 


1.0 


2.447376 


n/a 


1.5 


1.748874 


n/a 


2.0 


1.430883 


n/a 


2.6 


1.238051 


1.238012(#) 


3.0 


1.164321 


1.164319 


3.5 


1.105220 


1.105196 


4.0 


1.068256 


1.068237 


5.0 


1.029356 


1.029348 



Table 7: The first excited state of the scaling Lee- Yang model. The energy and the volume 
are normalized to the mass of the lowest excitation, which is the first breather of the unrestricted 
sine-Gordon model. The TCS data shown have the predicted bulk energy term subtracted. 



5.1 The Vir(2,2n+1) + $(1,3) series 

Let us start with examining the scaling Lee- Yang model Vir(2,5) + $(1,3). There is only 
one independent value of the twist which we choose as 

71 

since we have a single ground state, and as a result there are no kinks in the spectrum. 
We fix the value of a as above, so we still have a freedom of choosing 5. This can be done 
by matching to the UV dimensions: if for a certain state we choose the wrong value of 5, 
we find a conformal dimension that is not present in the Kac table of the model. 

The excited states are multi-particle states of the first breather of the corresponding 
unrestricted sine-Gordon model, which has 

2 

Now one can calculate the state containing one particle at rest. We find the numerical 
data presented in table |?]. 

It turns out that as we decrease I, the self-conjugate root starts moving to the right. 
It does not remain in the middle like in the a = case, which is to be expected since for 
nonzero a we have no left/right symmetry. However, the total momentum of the state still 
remains zero due to a contribution from the integral term in eqn. ( |2.12|) . One can see that 



once again we have the phenomenon noticed in the case of the first breather of sine-Gordon 
theory, namely the appearance of the special root and its two accompanying holes, so the 
iteration breaks down again around I = 2.5. The (#) in the table |7] written after the NLIE 
result for I = 2.6 means that due to the fact that the singularities corresponding to the 
new holes and the special root are just about to cross the contour and upset the iteration 
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scheme, the NLIE result becomes less precise. We will use this notation on later occasions 
too. In any case, the agreement still looks quite convincing. 

Let us now look at the UV spectrum of the model. We know that the Lee- Yang model 
contains only two primary fields, the identity I and the field (p with left/right conformal 
weights 

A+ = A" = -- . 

5 

In fact, the ground state of the massive model corresponds to (p in the UV limit. One 
can compute the UV limit of the first particle from the NLIE too, taking into account the 
appearance of the special root and the holes. It turns out that the special root and one 
of the holes moves to the left together with the self-conjugate root, while the other hole 
moves to the right. The result is 

A+ = A" = , 

i.e. the identity operator I, which fits nicely with the TCS data (see also |T6| ). 

Let us look now at moving breathers. If the self-conjugate root has Bethe quantum 
number 1=1, the corresponding state will have momentum quantum number 1, i.e. 

P = ^ 
R ' 

and in the UV A + — A~ = 1. One can note from the numerical data presented in table 
H that the special root does not appear here. The reason is that the self-conjugate root 
moves to the left and the real part of its position -d is given to leading order by 

2ttI 

sinh(9fte$) ~ — - — . 

Ib 

As a result, the contribution to the derivative of Z from the Zsinhi? term remains finite 
when I — ► 0. In the previous example of the particle at rest the left-moving nature of the 
self-conjugate root when 1 = does not prevent the occurrence of the breakdown in the 
iteration scheme: since its Bethe quantum number is zero, it does not move fast enough 
to the left in order to balance the negative contribution to derivative of Z coming from 
the self-conjugate root source. At the moment we have no way of predicting analytically 
whether or not there will be specials in the UV limit: we just use the numerical results to 
establish the configuration for the evaluation of UV weights, supplemented with a study 
of the self-consistency of the solution of the plateau equation ( |L33| ). The UV dimensions 



for the moving breather turn out to correspond to the state L_xip. 

One can similarly compute the UV dimensions for some other excited states. For 
example, the two-particle states with half-integer Bethe quantum numbers I\ > 0, I2 < 
for the two self-conjugate roots are found to have 

A+ = -- + h + - , A- = -- - h + - , 

5 2' 5 2 ' 
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Ib 


TCS 


NLIE 


0.1 


60.75048 


60.74682 


0.5 


12.20618 


12.20561 


1.0 


6.182516 


6.182363 


1.5 


4.202938 


4.202915 


2.0 


3.231734 


3.231640 


2.5 


2.662186 


2.662110 


3.0 


2.292273 


2.292231 


3.5 


2.035552 


2.035530 


4.0 


1.848892 


1.848849 


5.0 


1.599792 


1.599762 



Table 8: The one-particle states with Lorentz spin 1 of the scaling Lee- Yang model. The energy 
and the volume are normalized to the mass of the lowest excitation, which is the first breather of 
the unrestricted sine-Gordon model. The TCS data shown have the predicted bulk energy term 
subtracted. 



1 


TCS 


NLIE 


0.1 


123.583 


123.5693 


0.5 


24.7806 


24.77936 


1.0 


12.4870 


12.48635 


1.5 


8.42926 


8.428693 


2.0 


6.42931 


6.429201 


3.0 


4.48444 


4.484209 


4.0 


3.56367 


3.563519 


5.0 


3.04899 


3.048881 



Table 9: The lowest lying zero-momentum two-particle state in the scaling Lee- Yang model as 
computed from the NLIE and compared with TCS. 



in agreement with TCS data which show that they correspond in the UV to descendent 
states of ip. The first such state with quantum numbers 

1 2 ' 2 

corresponds in the UV to L_iL_i<£> and is given numerically in table ||. 

The lowest lying three-particle state of zero momentum, with Bethe quantum numbers 
(—1,0, 1) corresponds to the left/right symmetric second descendent of the identity field, 
i.e. to the field TT, where T denotes the energy-momentum tensor. This is very interesting, 
since from experience with NLIE UV calculations one would naively expect this to be a 
first descendent (descendent numbers are usually linked to the sum of Bethe quantum 
numbers of left/right moving particles - see the examples in -, and this state is the 



lowest possible descendent of the identity I). However, the field L_\L_xL is well-known to 
be a null field in any conformal field theory. 
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The above correspondences are again confirmed by comparing to TCS (see the wonderful 
figures in Ql6j| ) . In general, one can establish the rule that states with odd number of 
particles must be quantized by integers (5 = 1), while those containing even number of 
particles must be quantized by half-integers (5 = 0) in order to reproduce correctly the 
spectrum of the scaling Lee- Yang model. 

We conducted similar studies for the models Vir(2, 7) + $(1,3) and Vir(2,9) + $(1,3) 
and found similarly good agreement with TCS data. For the first one-particle state of the 
model Vir(2, 7) + $(1,3) we also checked our results against the TBA data in the numerical 
tables of |3l| and found agreement with the TBA results. 

Given the choice of a above, the correct rule of quantization in all of the models 
Vir(2,2n + 1) + $ ( i )3) is 

5 = M sc mod 2 , 

where M sc is the number of self-conjugate roots in the source corresponding to the state. 



This is exactly the same rule as the one established for pure sine-Gordon theory in p5| . In 
the presence of the twist, such a rule of course has meaning only together with a definite 
convention for the choice of a. 



5.2 One-breather states in the Vir(3, 7) case 

It is interesting to note that in the case of Vir(3,n) models, all the neutral states must 
come in two copies, since they can be built on top of either of the two ground states. We 
take the example of the Vir(3, 7) model and the states corresponding to a breather at rest. 
We have 

3 

but now there are two inequivalent values for the twist 

7T 2n 
01 ~ 3 ' T' 



When a = 0, we can calculate the critical value of I to be l cr iucai = 5.23 using (|2.19| ). In 



this case, the twist helps a bit, because it makes the self-conjugate root a left mover; it is 
intuitively clear that the bigger the twist, the more it lowers the eventual value of l cr ucah 
which is in accord with the numerical results of table [R]. From the TCS data one can 
identify that breather #1 is really the one-particle state in the sector of ground state #1, 
while breather #2 is in the sector over ground state #2 using the notations of table ||. 
A direct calculation of the conformal weights gives the following results: 

A+ = A- = — 
28 

for breather #1 and 

A + = A~ = 
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breather 1 


breather 2 


1 


TCS 


NLIE 


TCS 


NLIE 


0.1 


32.21645 


n/a 


18.76030 


n/a 


0.5 


6.671964 


n/a 


4.037716 


n/a 


1.0 


3.662027 


n/a 


2.434501 


2.434431(#) 


1.5 


2.769403 


n/a 


2.027227 


2.027213 


2.0 


2.385451 


n/a 


1.884404 


1.884388 


2.5 


2.190232 


n/a 


1.829459 


1.829456 


3.0 


2.079959 


n/a 


1.809244 


1.809248 


3.5 


2.012596 


n/a 


1.803873 


1.803886 


4.0 


1.968808 


1.968784(#) 


1.804937 


1.804953 


4.5 


1.938895 


1.938889 


1.808633 


1.808658 


5.0 


1.917648 


1.917635 


1.813191 


1.813224 



Table 10: The two one-breather states of the Virasoro minimal model Vir(3,7) perturbed by 
$(i,3)- The energy and the volume are normalized to the mass of the kink, which is the soliton of 
the unrestricted sine-Gordon model. Breather #1 has a = |, while breather #2 corresponds to 
a = The TCS data shown have the predicted bulk energy term subtracted. 



for breather #2, which are in complete agreement with the TCS data. We also checked 
the two different states containing two breathers with Bethe quantum numbers 




and found an equally excellent numerical agreement with TCS. Just like in the case of 
sine-Gordon and scaling Lee- Yang model, for these states one can continue the iteration 
of the NLIE down to any small value of I, although at the expense of a growing number of 
necessary iterations to achieve the prescribed precision. 



6 Conclusions 

In this paper we described how to use the NLIE approach to compute finite size effects in 
Virasoro minimal models perturbed by $(1,3). Up to now, apart from a previous attempt 
by P. Zinn- Justin Hl9| , the only results in this framework were the description of ground 



states of $(1,3) perturbed unitary minimal models |TT . 

As a starting point, we first extended our previous studies |14|, [15], [17| of sine-Gordon 
theory to breather states which then played a prominent role in providing us concrete 
examples of excited states in minimal models. We showed that in the IR limit the NLIE 
successfully reproduces the scattering amplitudes involving breathers and we gave examples 
of comparing numerical results from the NLIE to those coming from TCS. 

We then proceeded to the case of perturbed minimal models. Using a twisted NLIE a la 
Zamolodchikov [12|] we provided examples for ground states in nonunitary models, thereby 
overcoming the limitation of flTl| which considered only the unitary case. 



29 



Concerning excited states, we started by establishing the connection of the twisted NLIE 
to a twisted version of the lattice Bethe Ansatz. Doing so, we have found two shortcomings 
of the results in The first of these was an incorrect relation (|4.28| ) between the twist 



parameter u in the lattice Bethe Ansatz and the one appearing in the NLIE, which did not 
reflect the periodicity of the Bethe Ansatz in the u, in contrast to the relation ( |4.27| ) found 



by us. The other was that to describe all possible states of the UV limiting minimal CFT, 
it is not enough to choose just one value of the twist parameter: indeed a whole range of 
values is required ( |3.23| ). We demonstrated that the twisted NLIE gives conformal weights 



that are consistent with the spectrum of minimal models. Numerical calculation of concrete 
examples gave us a strong evidence for the correctness of the energy levels derived from 
the twisted NLIE for excited states. 

Throughout this paper we made one important omission: we did not treat the multi- 
kink states so characteristic of the perturbed minimal models. Although the general treat- 
ment of the UV spectrum in section 4.2 is valid for those states too, a detailed description 
is far more complicated than for states which contain only breathers and is left open to 
further studies. 

We also pointed out a technical difficulty, namely that the source configuration of the 
NLIE may change as we vary the volume parameter I. Typically what happens is that while 
the counting function Z is monotonic on the real axis for large volume, this may change as 
we lower the value of / and so-called special sources (and accompanying holes) may appear. 
We do not as yet have any consistent and tractable numerical iteration scheme to handle 
this situation, although the analytic UV calculations and intuitive arguments show that 
the appearance of these terms in the NLIE is consistent with all expectations coming from 
the known properties of perturbed CFT. In addition, in the range of I where we can iterate 
the NLIE without difficulty, our numerical results show perfect agreement with TCS. We 
want to emphasize that these transitions are not physical: the counting function Z and 
the energy of the state is expected to vary analytically with the volume, it is just described 



by a NLIE with modified source terms. As it was pointed out in [|T^, |ToJ, the whole issue 
is related to the choice of the branch of the logarithmic term in the NLIE ( ]2.5|) . The 
problem itself is very similar to the behaviour of singularities encountered in the study of 
the analytic continuation of the TBA equation [|4j and we can hope that establishing a 
closer link between the two approaches can help to clarify the situation. 

We think that the work presented in this paper provides a strong evidence for the NLIE 
description of excited states in perturbed minimal models. We would like to mention that 
following the lines of [H]] it is possible to extend this framework to minimal models of 
W algebras based on the ADE Lie algebras. It seems very likely that apart from the 
loopholes (concerning the twist parameters) pointed out and clarified in this paper no 
further complications will occur, but a detailed discussion of these issues is out of the 
scope of the present publication. 
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